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Abstract 



■ A is a set of N vectors in Z N ~ 2 situated in a hyperplane not through 

and spanning Z N ~ 2 over Z. Gulotta's algorithm [3] constructs from A 
a dimer model. A theorem in [6] states that the principal A-determinant 
• equals the determinant of (a suitable form of) the Kasteleyn matrix of 

that dimer model. In the present note we translate Gulotta's pictorial de- 
scription of the algorithm into matrix operations. As a result one obtains 
an algorithm for computing the principal A-determinant, which is much 
faster than the algorithm in [5]. 



1 Introduction 

00 - 

\Q . A = {ai, . . . , ajv} is a set of vectors in Z situated in a hyperplane not 

fT") ' through and spanning Z N ~ 2 over Z. The principal A-determinant, de- 

fined in [2], describes the singularities of Gelfand-Kapranov-Zelevinsky's 
^ ■ A-hypergeometric system of partial differential equations [3]- It also de- 

' scribes for which iV-tuples of coefficients u\ , . . . , ujv £ C the Laurent 

polynomial X^^Li u i xaj m N — 2 variables is singular (see [2] for details). 
It is a polynomial with integer coefficients in the variables ui, . . . ,un- 
The restriction rankA = $A — 2 means that the corresponding hypergeo- 
, metric functions are essentially functions in two variables and that in the 

5^ ' Laurent polynomial the number of terms exceeds the number of variables 

by 2. Even in this case the definition of the principal A-determinant is 
fairly complicated and only a few of its coefficients could explicitly be 
calculated in [2]. In [T] Dickenstein and Sturmfels re-examined the defini- 
tion of the principal A-determinant and related it to Chow forms which is 
another important concept from [2]. In [6] it was shown how these Chow 
forms and the principal A-determinant can be easily computed as the de- 
terminant of a suitable version of the Kasteleyn matrix of a dimer model 
associated with A. When writing 6 I had only the algorithm in [5] to 
construct that dimer model. Later Gulotta gave another algorithm [4] for 
constructing the appropriate dimer model. He describes the algorithm as 
a process that transforms certain doubly periodic configurations of curves 
in the plane. In Sections [2] [3] |4] we give a faithfull reproduction of those 



1 



configurations of curves by matrices and of Gulotta's algorithm by row 
and column operations on these matrices. In that form Gulotta's algo- 
rithm, which is a fast converging iterative process, is much more efficient 
than the algorithm in [5], which is a search with many trial-and-errors. 
Moreover, unlike for the algorithm in [5] for Gulotta's algorithm it can be 
guaranteed that it finds a desired dimer model. On the other hand there 
are cases in which [5] yields two different models and [4] gives only one. 

In this note we do not need formal definitions of 'dimer model' and 
'principal ^-determinant'. Dimer models are implicitly present through 
their patterns of zigzags. This is briefly explained in Remark 12.41 Princi- 
pal ,4-determinants appear only in Section |S] in a quotation from [B]. 

In Section [5] we recall from [5] how a pattern of zigzags (alias dimer 
model) is faithfully represented by a matrix Kz(z,u), which is in fact a 
suitable generalization of the Kasteleyn matrix of the dimer model. In 
Section [6] we recall from [6] the theorem that expresses the principal A- 
determinant as the determinant of K,z(z, u). Sections [5] and [6] and can be 
read immediately after Definitions 12.21 and 12.31 

2 Patterns of zigzags on the torus 

2.1. This note is about patterns of zigzags on the torus T = R 2 /Z 2 . Here 
zigzag means the image (modulo Z 2 ) of an oriented connected curve C in 
the plane (i.e. the image of a continuous map from R to R 2 ) such that the 
two coordinate functions restrict to monotone functions on C and such 
that t + C — C for some non-zero vector t £ Z 2 . 

We denote by e\ (resp. £2) the zigzags in T coming from the first 
(resp. second) coordinate axis of R 2 . The homology classes of ei and 
£2 form the standard basis for the homology group i?i(T, Z). We denote 
the intersection number of two elements a, f3 £ H\(T,Z) by a A /3. For 
two zigzags Z, Z' on T we write Z A Z' for the intersection number of 
their homology classes. For two zigzags Z, Z' which intersect in only a 
finite number of points and for which the intersections are transverse, 
each intersection point contributes +1 or —1 to Z A Z' according to the 
orientation; e.g. £1 A £1 = £2 A £2 = 0, £1 A £2 = — £2 A £1 = 1. 

Writing [Z] for the homology class of a zigzag Z we have 

[Z] = (ZAEa)[ei] - (ZAei)[e 2 ]. (1) 

2.2. Definition. In this note pattern of zigzags means a finite sequence 
Z — (Zi, . . . , Z p ) of zigzags on T which satisfies the following conditions: 

1. The homology classes [Z{\, . . . , [Z p ] span Hi(T, Q). Their sum is 0. 

2. For every j the greatest common divisor of (Zj A £1, Zj A £2) is 1. 

3. Every point of T lies on at most two zigzags in Z. 

4. Every pair of zigzags Zi, Zj intersects in only a finite number of 
points and the intersections are transverse. 

5. The 2-cells (i.e. connected components) of T \ Uf=i are divided 
into three types: those (called +- cells) of which the boundary is 
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positively oriented, those (called — cells) of which the boundary is 
negatively oriented and those of which the boundary is not oriented. 
It is required that for every intersection point x of zigzags in Z one 
of the four 2-cells having x in their boundary is a +-cell, one is a 
— cell and two have an unoriented boundary; see Figures 111 1UI 

6. The number of +-cells equals the number of — cells. 

2.3. Definition. We say that a pattern of zigzags Z — (Zi, • ■ • , ^v) is 
good if in addition to the above conditions, it also satisfies: 

7. Write Q for the column vector [Zi A ei, Zi A £2] • Then the de- 
terminants det(£i,Ci+i) for i — 1, ... ,p — 1 and det(£ p ,Ci) are not 
negative. Moreover, Q 7^ — Ci+i for i = 1, . . . ,p — 1 and £ p 7^ ±£1. 

8. Every ordered pair (Zi, Zj) of zigzags has the same orientation at 
all its intersection points. This is equivalent with: 

Vz,j : \Zi AZj\ = KZ.nZj). 

Condition 7 means that the homology classes are ordered counterclock- 
wise with increasing indices and that for every homology class the indices 
of the zigzags in that class form a connected interval in the index set 
{1, ... ,p}. In l2.12l we formulate a condition which also restrains the or- 
dering of zigzags within their homology class. 

2.4. Remark. The dimer model for a pattern of zigzags is the graph 
with a node • for every +- cell and a node o for every — cell. Two nodes 
are connected by an edge if the cells have a vertex in common. Figure [T] 
shows an example. 




Figure 1: Zigzag pattern and corresponding dimer model. 
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2.5. To a pattern of zigzags Z = (Zi, 



, Z p ) we assign the 2 x p-matrix 



Z\ A ei , . . . , Z p A Ei 
Zi A e 2 , ■ • ■ , Z p Ae 2 \ ' 

displaying the intersection numbers of the zigzags in the pattern with the 
two curves ei and £2- We assume that all intersections of zigzags with £1 
and £2 are transverse. The intersection numbers are visible in the pictures 
as follows. Represent T by the unit square with opposite sides identified. 
Then Zj A £2 is the number of times the zigzag Zj crosses the right-hand 
vertical edge from left to right minus the number of times it crosses from 
right to left. Similarly, Zj A £1 is the number of times Zj crosses the top 
horizontal edge downwards minus the number of times it crosses upwards. 
Condition 12.21 1. is equivalent with 

the rank of Bz is 2 and the sum of its columns is 0. 

2.6. From fTJ we see: 

ZiAZj = (ZiAei)(ZjA£2)-(ZiAe 2 )(ZjAei) = det 

This together with Condition 12.21 2. implies 

ZiAZj =0 [Zi] = ±[Z S ]. 

And thus, rank.Bz = 2 if and only if [Zi] 7^ ±[Zj] for some i,j. 

2.7. Let Z — {Zi, . . . , Zp} be a pattern of zigzags on T. Pick a point * 
in one of the — cells. To every 2-cell c we associate a row vector in V as 
follows. Take any path 7 on T starting at * and ending in (the interior 
of) c, such that 7 intersects zigzags transversely. Each point in Zj n 7 
contributes, depending on the orientation, +1 or —1 to the intersection 
number Zj A 7. Then to c we associate the vector of intersection numbers, 
or briefly intersection vector, [Z\ A 7, . . . , Zp A 7]. Choosing another path 
7' from * to c changes this vector by a Z-linear combination of the rows 
of the matrix Bz- 

It follows from Condition l2.1l 5 that at an intersection point of zigzags 
Zi and Zj the vectors for the two cells with unoriented boundary and the 
cell with negatively oriented boundary can be obtained from the vector for 
the +- cell by subtracting 1 from the i-th coordinate, respectively 1 from 
the j-ih coordinate, respectively 1 from both 2-th and j-th coordinate. We 
can thus capture all relevant information of the pattern of zigzags Z in 
the matrix Bz and two additional matrices Iz and Pz, denned as follows. 

2.8. Definition. The columns of Iz and Pz correspond with the zigzags 
Z\, . . . , Zp. The rows of Iz and Pz correspond with the intersection points 
of pairs of zigzags in Z. The row of matrix Iz for a point x S Zi PI Zj 
has 1 in positions i and j and elsewhere. Matrix Pz has in the row for 
intersection point x an intersection vector of the +- cell which has x in its 
boundary; see Figures [21 and ITUl for examples. 



B z = 



(ZiA£i) (ZjAei) 
(Zi A £2) (Zj A £2) 
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It is also convenient to have the short notation Qz = Pz — Iz- Then 
matrix Q z has in the row for an intersection point x an intersection vector 
of the — cell which has x in its boundary. 
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Figure 2: Pattern of zigzags and the corresponding matrices. 



2.9. Remark. Due to the Z 2 i?2:-ambiguity in the choice of the intersec- 
tion vectors there is also a Z 2 _B^-ambiguity in the rows of the matrices 
Pz and Qz in Definition 12,81 In the algorithm we start with well-defined 
matrices Pz and Qz- In the course of the algorithm we only delete rows 
and perform the same operations on the columns of the matrices Bz, Pz, 
Qz and Iz simultaneously. So the algorithm is also unambiguous. It does 
however happen that rows of Pz (resp. Qz) which corresond to the same 
+- cell (resp. — cell) are not the same, but differ by a vector in I?Bz- 

In Section [5] we pass to V/l?Bz- 

2.10. Notation. For a matrix M we denote j-th column as M(:,j), the 
i-th row as M(i, :) and the (i, j)-entry as M(i,j). 

2.11. Definition. Let Zj and Zk be two zigzags such that [Zj] = — [Zk]- 
We say that (Zj, Zk) is a -^-opposite pair (resp. —-opposite pair) if 

Qz(:,j) = -Qz(:,k) (resp. P z (:,j) = -Pz(:,k)). 

Suppose Zj (~l Zk = 0. Then (Zj, Zk) is a +-opposite pair (resp. — 
opposite pair) if and only if there are between Zj and Zk no — cells (resp. 
no +-cells). Most pictures in this note contain examples of opposite pairs. 
The term 'opposite pair' without ± was introduced in [4 §5.3. 

2.12. Definition. A good pattern of zigzags Z — (Z\, . . . , Z v ) is said to 
be very good if it satisfies: 

9. For every homology class [Z] the sequence of all zigzags (Zi, . . . , Zi+ r ) 
in homology class [Z] and the sequence of all zigzags (Zj, . . . , Zj+ S ) 
in homology class —\Z\ satisfy: 
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for < t < min(r, s) : 
and either: 

for < t < min(r, s) 
or: for < t < min(r, s) 

3 The moves in the algorithm 

Gulotta's algorithm transforms in an iterative way a very good pattern of 
zigzags Z into another very good one Z' . In [4] the algorithm is mainly 
described by transforming a drawing of Z into a drawing of Z' . We will 
present the same algorithm by row and column operations on the matrices 
B z , Iz, Pz- 

3.1. Merging move. The basic move in the algorithm merges two zigzags 
Zi and Zj which intersect in exactly one point, as shown in Figure^ 



(Zi+t, Zj+t) is a -{--opposite pair 

(2) 

(Zi+t+i, Zj+t) is a —-opposite pair 
: (Zi+t, Zj+t+i) is a — opposite pair 



Figure 3: Merging move. 



It is evident that the merging moves preserve Conditions 1-6 in !2,2l 

It was pointed out in |4], that when two zigzags Zi and Zj of a pattern 
Z are merged and become one zigzag Z in Z' then 



Z A £i 
Z A £2 



ZiAei 
ZiAe 2 



+ 



Z 3 A £i 

Zj A £2 



(3) 



Actually this means [Z] = [Zi] + [Zj] and, hence, also the column of the 
matrix P z i which corresponds with the zigzag Z is the sum of the i-th 
and the j'-th columns of Pz- Moreover, as the picture indicates, the point 
of intersection Zi n Zj disappears. The following statement also specifies 
where we put the new zigzag Z in the list of zigzags for Z' . 
Conclusion: The merging of Zi and Zj for Zi A Zj — 1 is given by the 



same column operation on Bz, Iz, Pz, namely: 
the i-th and subsequently delete the j-th column, 
and Pz the row for Zi n Zj. 



add the j-th column to 
It also deletes from Iz 



Merging moves performed on a very good pattern of zigzags need not 
preserve Conditions 12.31 7-8 and 12.121 9. Some repairing may be needed in 
order to turn the pattern of zigzags produced by the merging moves into 
a very good one again. 



3.2. Repairing move 1. The first type of repairing move is shown in 



Figure^ This is used when [Zi] — 
area between the zigzags is just one 



■[Zj] and nZj) = 2, while the 
—cell (as suggested in the picture) 
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Figure 4: Repairing move 1. 



or one +-cell (interchange + and — in the picture) . From this picture one 
immediately comes to the conclusion: 

Conclusion: Repairing move 1 just deletes the rows for the two points of 
Zi n Zj from Iz and Pz ■ 

3.3. Repairing move 2. The second type of repairing move is shown 
in Figure^ This is used when [Zi] = [Zj] and %{Zi PI Zj) = 2. In this 




Figure 5: Repairing move 2. 

case one may distinguish three kinds of rows in the matrix Pz, according 
to whether the i-th entry minus the j-th entry is equal to m — 2, m — 1 
or m, for some integer m (depending on i and j). The rows of the latter 
kind correspond in the picture with +-cells in the area between the two 
zigzags. When the intervals of Zi and Zj between the points of Zi n Zj 
are swopped (as suggested by the right-hand picture) one must add 1 to 
the j-th coordinate and subtract 1 from the i-th coordinate in all rows 
of Pz corresponding with a +- cell in the area between Zi and Zj. One 
must also interchange the i-th and j-th entries in the rows of Iz which 
correspond with intersections with the intervals of Zi and Zj between the 
two points of Zi l~l Zj . And one must delete from Iz and Pz the two rows 
corresponding to the two points of Zi n Zj . 

Conclusion: Repairing move 2 operates on the columns of Pz as follows: 
Write H(r) = Pz(r,i) — Pz(r,j) and m = max r (_H"(r)). Then 

Pz(r,i) ~* P z {r,i)-1, Pz{r,j) P z (r,j) + 1 if H(r)=m, 
Pz(r,i) - Pz(r,i), Pz(r,j) Pz(r,j) if H(r)^m. 



7 



It operates on the columns of Iz by: 

Iz(r,i) •v* Iz(r,j), Iz(r,j) ~» Iz(r,i) 

if H(r) — m and Iz(r,i) = 1, or H(r) =m — l and Iz(r,j) = 1, 
Iz(r,i) ~-» Iz(r,i) , Iz(r,j) -*+ Iz(r,j) otherwise. 

Finally, it deletes from Iz and Pz the rows for the points of Zi n Zj . 



12 3 4 1" 2" 3" 4' 




Figure 6: Transforming an alternating sequence of opposite pairs. 

3.4. Example. Let {Z\, . . . , Z% 3 ) and (Z[, . . . , Z^ s ) be two sequences 
of zigzags in a pattern Z' 1 ' such that %{Zi (~l Zj) = 1 for all i,j and 
Zi n Zj = Zl C\ Zj = $ for all i / j and such that (Zj-i,Zj) and 
(Z / J _ 1 , Zj) are ( — l^-opposite pairs for j — 2, ... ,2s. There are two well 
controllable cases in which merging of these two sequences followed by 
repairing moves 2 and 1 yields a sequence of zigzags (Z'l" , . . . , Z'^'J) such 
that Z'l" n Z'J" = for all i / j and such that {Z'!'\,Z$") is a 
opposite pair for j = 2, . . . , 2s. 

Case 1: Zj and Zj merge for j = 1, ... , 2s. 
Case 2: and Zj_/ 1 p merge for j = 1, . . . , 2s. 

Figures [6] and [7] show this for s = 2 and clearly generalize to arbitrary s. 

In either case let (Z'l , . . . , Z'l s ) be the sequence of zigzags in the pattern 
ij' 2 ' which results from the merging. Now transform Z^ by applying 
repairing moves of type 2 at the points of Z" n Zj for all i ^ j for which 
[Zl'] = [Zj]. The result is the sequence of zigzags (Z[" , . . . , Z'l' s ) in the 
pattern Z {S) . Then Z'l' n Z"' ^ only if [Zf] = -[Zj"]. Next transform 
.Z' 3 ' by applying repairing moves of type 1 at the points of Z"' HZ"' for all 
i, j. Call the resulting pattern of zigzags Z^ and the relevant sequence of 
zigzags (Z'{", Z'l"). In this last sequence Z'l" n Z"" = for all i / j 
and (Z'j\,Z'j") is a (-l) 3 -opposite pair for j = 2, . . . , 2s. 
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1 2 3 4 1" 2" 3" 4' 




Figure 7: Transforming an alternating sequence of opposite pairs. 

It is instructive to perform the moves in Figures [6] and [7] also for the 
matrices Bz, Iz, Pz- For the top-left picture Bz, Iz, Pz are given in 
Figure 1101 For the other pictures one may follow the description of the 
merging and repairing moves. 

3.5. Remark. In 13.41 we have chosen the labels for the zigzags while 
drawing the pictures. If one uses the matrix operations instead, the al- 
gorithm determines the labels and it may be necessary to reorder the 
columns of the matrices Pz, Qz, Iz to meet the requirements of Equa- 
tion @. Example 13.41 shows that such a reordering is always possible. 

3.6. Repairing move 3. The third type of repairing move is shown in 
Figure^ It is used for a zigzag Zo and a sequence of zigzags (Zi, . . . , Z2s) 
which satisfy the following conditions. Firstly, ZiClZj — for all z > j > 1 
and (Zj-i, Zj) is a ( — l) J -opposite pair for j = 2, ... , 2s. Secondly [Zo] = 
i[Z\] and (j(Zo n Zf) — 2 for j = 1, . . . , 2s. Reversing if necessary the 
labeling in the sequence (Zi, . . . , Z2s) we may without loss of generality 
assume [Z ] = (-l) 3 ^] for j = 1, . . . , 2s. 

Gulotta's instructions (cf. [4] §5.3) in this situation are to remove 
(Zi, . . . , Z2s) and to insert a sequence of zigzags (Z[ , . . . , Z' 2s ) such that 
Z[ nZ'j = for all j > i > and such that (Z' j _ 1 , Zj) is a ( — l)-?-opposite 
pair for j = 1, . . . , 2s. For notational convenience we write here and below 
Z' = Zo. 

In terms of the matrices Iz and Pz this means that we first delete 
from Iz and Pz all rows which correspond with an intersection point on 
one of the zigzags Z\, . . . , Z 2s and subsequently replace, for j = 1, . . . , 2s, 
the column of Pz which corresponds with the zigzag Zj by ( — 1) J times 
the column of Pz which corresponds with the zigzag Zo- 

Next we expand every row of Iz and Pz which corresponds with an 



9 



z 




Figure 8: Repairing move 3. 



intersection point of Zo and a zigzag Z x ^ Zo, Z\, ■ ■ ■ , Z2s to 1 + 2s rows 
which correspond with the intersection points of Z^ with Z' , Z[, . . . , Z' 2li 
(see Figure [9} . The columns of Iz and Pz are labeled in such a way that 
the column which originally corresponded to Zj now corresponds to Z'j 
for j = 0, . . . , 2s. 



r?t rjl ry! ry / ry I 

^0 ^1 ^2 ^2s-1^2s 



Figure 9: Intersections with alternating sequence of opposite pairs. 



Thus a row of Iz for an intersection point x of Z^ and Zq is replaced 
by 1 + 2s rows of which the j'-th one (for j = 0, . . . , 2s) has entry 1 in the 
columns corresponding with the zigzags Z^ and Z'j. Of the 1 + 2s new 
rows of Pz the 0-th one is equal to the row of Pz which corresponds to x. 
Figure [9] shows that for odd j the j-ih row is obtained from the (j — l)-st 
one by adding —1 in the column for Z'j_ x and +1 in the column for Z'j\ 
for even j > 2 the j-ih row is equal to the (j — l)-st one. 



4 Running the algorithm 

In this Section we translate the algorithm described by Gulotta in terms 
of pictures, into an iterative proces operating on matrices. 
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4.1. In order to prepare the input for the algorithm from the set A = 
{ai, . . . , ajv} (see JTJ we take a 2 x N- matrix Ba such that its rows are a 
Z-basis for the lattice {(h, ■ ■ ■ , In) € 1 N \ li&i + . . . + In&n = 0}. 

The aim of the algorithm is to create a very good pattern of zigzags Z 
such that Bz results by permuting and splitting up the columns of Ba as 
follows 



column 



of Ba ~+ d — g.c.d.(r, s) columns i 



in B z - (4) 



4.2. The algorithm starts with a very good pattern of zigzags Z for which 
Bz is the 2 x (2m + 2ri2)-matrix 



1...1 0...0 -1 1 0... 

0...0 1...1 0... -1 1 

where n\ (resp. 712) is the sum of the positive entries in the first (resp. 
second) row of Ba- This matrix is realized by a pattern with straight lines, 
m vertically down, m vertically up, 712 horizontally left-to-right and ri2 
horizontally right-to- left. To get a very good pattern one takes the vertical 
lines alternatingly down and up, and the horizontal lines alternatingly 
left-to-right and right-to- left. The matrices Bz, Pz and Iz for the initial 
pattern have 2ni + 2772 columns. There are 4niri2 intersection points and, 
hence, Iz and Pz have 4ni7i2 rows. We build Iz and Pz as follows: the 
+- cells are given by pairs (a, b) in {1, . . . , n\} x {1, . . . , 712}. The rows of 
Iz which correspond to the four vertices of the +- cell (a, b) are have 1 in 
positions a, b + ni, resp. a, b + 2ni + 71,2, resp. a + n\ + 712, & + m, resp. 
a + ni+77,2, b + 2ni + 712 . All other entries in these rows of Iz are 0. The 
non-zero entries of the four rows of Pz which correspond to the vertices 
of the +- cell (a, b) are 1 in position j 'if 1 < j < a or 2n\ + n2 + 1 < j < 
2rti + 772 + b and —1 in position j if m + 712 + 1 < j < ni + n>2 + a — 1 or 
ni + 1 < j < m + b — 1. Figure [lOl shows an example with ni = 772 = 2. 

4.3. At the beginning of an iteration step we have a very good pattern of 
zigzags Z — {Z\, . . . , Z p }. The merging moves are determined from the 
positions of the columns of the matrix Ba with respect to the columns of 
the matrix Bz- The columns of Bz are vectors in the plane R 2 and the 
ordering by increasing index coincides with the counter-clockwise cyclic 
ordering. In agreement with this cyclic structure we treat the first column 
of Bz as consecutive to the last column. 

A column v of Ba is either a positive integer multiple of a column of 
Bz or there is a unique pair of consecutive columns vi and 7J2 of Bz such 
that v = C1W1 + C2W2 with ci, C2 £ Q>o- 

4.4. The algorithm terminates automatically when all columns of 
Ba are multiples of columns of Bz- Since the merging moves decrease 
the number of zigzags the algorithm will surely terminate. 
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Figure 10: Pattern to start from. 



4.5. Cramer's rule explicates the relation v — c\v\ + C2V2: 

det(ui, 112) u = det(i>, V2) V\ — det(i>, vi) V2 . (5) 

In the pattern we start with the determinants of consecutive pairs of non- 
equal columns of Bz are 1. The merging of two zigzags in a very good 
pattern Z with exactly one intersection point replaces the corresponding 
columns of Bz by their sum. Thus in the next iteration step in the 
algorithm Equation ([5]) becomes 

det(t>i, V2) v — (det(«, 172) — m) «i + m(vi + V2) — (det(v, v\) + m) V2 

with m = min(det(i>i, v), det(n, i>2))- This then gives v either as a multi- 
ple of v\ + V2 or as a positive linear combination of vi and v\ + V2 or of 
vi +V2 and V2- Note that Aet(v\,v\ +^2) = det(ui J rV2,V2) = Aet(v\,V2)- 
Conclusion: In all cases in which Equation is used in the algorithm 
det(vi,^2) = 1 and the equation actually reads 

v — det(v, V2) vi — det(v, vi) V2 . (6) 

Column v of Ba thus leads to the command that m zigzags of the pattern Z 
in the homology class corresponding with the column vi of Bz must merge 
with m zigzags in the homology class corresponding with the column V2; 
here m = min(det(«i, v), det(u, V2)). 

4.6. As Equation ([6]) indicates we need the determinants of the 2x2- 
matrices with first column from Ba and second column from Bz- These 
are simultaneously given as the entries of the N x p-matrix 

S = B A JBz with J = 

The columns of Ba correspond with the rows of S. One can implement the 
discussion in 14.31 and 14.51 for all columns of Ba simultaneously as follows. 
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Let S c = S(:, [2 : p, 1]) be the TV x p-matrix obtained from S by cyclically 
permuting the columns so that the first column comes in the last position. 
Let R be the N x p-matrix with (i, j)-entry 



Rio = + + ifS»i<0, 

R i3 = if Sij > . 

We define functions p and A on {1, ... , p} by 

N 

p(J) = E ^ . A 0') = - 1) if i > 1 , A(l) = p(p) . 

1=1 

Next we define for a homology class of zigzags [Z] in the pattern Z = 
(Zi, . .. , Z p ): 

p([Z}) = max{p(j)\Zj€[Z]}, 
\{[Z]) = max{A(j) | Zj G [Z]} , 
W[Z\) = «([Z]) - p([Z}) - ~X([Z}) . 

Then p([Z]) (resp. A([Z]) ) is the number of zigzags in [Z\ which must 
merge with a zigzag in the homology class immediately after (resp. before) 
[Z] and J1{[Z]) is the number of zigzags in [Z] which must not merge with 
another zigzag. The merging step defined in 14.81 decreases the number of 
zigzags by 

9 = X>(fc) 

h=l 

and uses the map <p : {1, . . . ,p} — > {1, . . . ,p — q}, 

m = { i-TLiMh). if i>A(i) (7) 

nj> I P-1 + J-J2UHh) if i<A(l). 10 



4.7. In order to eventually satisfy Requirement 12.121 9 and to benefit 
from Example 13.41 we permute the zigzags in each homology class as fol- 
lows. First we define for each homology class of zigzags [Z] for which the 
opposite class —[Z] also occurs in the pattern Z: 

p([Z]) = min{p([Z]),K-[Z])) , K[Z\) = mm{X([Z]),X(-[Z])) , 
M[Z]) = rmnmZ]),H-[Z])). 

If -\Z\ does not occur in Z we put p([Z]) = \{[Z]) = p.{[Z]) = 0. 
Next we write for every homology class [Z] in Z: 

p([z])=p([z])-p([z]), K[z])=K[z])-x([z]),H[z])=mz])-m])- 

The permutation we apply to the zigzags in homology class [Z] is a so- 
called shuffle. This means that \Z\ is split into disjoint intervals which 
are permuted, while inside each interval the ordering is unchanged. The 
shuffle we apply to the zigzags in [Z] is depicted in Figure [TT] the numbers 
A([Z]) etc. indicate the length of the interval. 
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x([z]) M[z]\ p([z}) Mz]l v{[Z\) , p([z]) 




\{[Z\) p([Z]) p([Z]) 

Figure 11: Shuffle within one homology class 



Such shuffles must be applied to each homology class in the pattern 
Z — (Zi, . . . , Z p ). The composite result is the shuffle permutation 

cr : {l,...,p} —»{!,..., p}. (8) 



4.8. Definition. We define the merging matrix Maz for the set A and 
the pattern of zigzags Z to be the p x (p — g)-matrix with (i, j')-entry 

{M AZ ) ij = lifj = ip{<r(i)), {Maz)h = if j / tp(v(i)) , (9) 

with <p and a as in © and ((8}. 

The merging step in the algorithm multiplies the matrices Bz, Iz, Pz 
and Qz from the right with the matrix Maz and subsequently deletes 
the rows which correspond to intersection points in the pattern Z which 
disappear in the merging process. These are recognized as the rows of 
Iz Maz with only one non-zero entry (namely 2). 

Thus the merging Z ~+ Z' is realized by 

' B z - = B z Maz 

I Z ' — delete rows from Iz Maz 
Pz> = delete rows from Pz Maz 
Qz- = delete rows from Q z Maz 

4.9. The transformation of patterns of zigzags Z Z' in (|10[) has 
been organized so that if Z is a very good pattern, then Z' satisfies the 
Conditions 12.21 1-6. 12.31 7 and the property formulated on the first line of 
Equation ([2]). However, Z' need not satisfy [2~3l 8. i.e. the equality 

IZlAZ'j] = %{Z[nZ'f). 

need not hold for all pairs of zigzags {Z[, Z'f) in Z 1 . This equality can only 
be violated if Z[ results from merging zigzags Zi t and Zi 2 from Z such 
that Z it A Zi 2 — 1 and (Z i± A Z'f)(Zi 2 A Z'f) < 0. Since [Z'j] can only be 
the homology class of a zigzag in Z or the sum of two such, this can only 
happen if [Z'j] = ±([Z h ] + [Z i2 ]) = ±[Z[]. 

The merging process in (|10|l is such that if Z[ results from merging 
zigzags Zi t and Zi 2 from Z, then every zigzag Z'j in the homology class 
[Z[\ (resp. in — [Z'i]) is the result of merging a zigzag Zj 1 from [Zi ± ] (resp. 
— [Zi x ]) with a zigzag Zj 2 from [Zi 2 ] (resp. —[Zi 2 ]). Since Z satisfies 
Condition 12.31 8. and Zi 1 A Zi 2 = 1 we have in that case 

z n n z n = z %2 n z ]2 = , \{z ix n z ]2 ) = %{z^ nz h ) = i. 
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Whence if Zj e ±\Z[\ and Z- / Zj, then %{Z[ n Zj) = 2. 
Conclusion: For a pair of zigzags (Z^Zj) in Z' we have: 

IZ-AZ'j] / %{Z[ n Zj) <^=> Z- A Zj = and %{Z[ n zj) = 2 . 

4.10. We now transform the pattern of zigzags Z' created in (|10|) into a 
very good one. First we apply repairing moves 2 fsee I3,3|l to those pairs 
of zigzags (Z|,Zj) which satisfy \Z[\ = [Zj] and \Z[ AZj| # {t(Z- n Zj) 
and which either both do or both do not belong to a +-opposite pair 
(cf. Definition ^. lip . Next we apply repairing moves 1 (see !3.2|) wherever 
possible. See also Example 1 3 . 4 1 and Remark [33] for the effect of repairing 
moves 2 and 1 on alternating sequences of opposite pairs. Thus with 
repairing moves 2 and 1 and possibly a permutation of columns in I z t, 
Pz>, Qz' we now have a pattern of zigzags which satisfies also the second 
half of Equation ((2J) and in which two zigzags Zi, Zj with [Zj] = ±[Zj] 
do not intersect unless one is member of an opposite pair and the other 
is not. 

Finally we apply repairing moves 3 (see !3.6p as follows for all homology 
classes [Z] and — [Z] which have resulted from merging and which satisfy 
tt([Z]) > [Z]) > 0. In these circumstances the zigzags in — [Z] and 
the first [Z]) zigzags in [Z] form a sequence (Zi, . . . , Zis) as in the 
beginning of 13.61 For the zigzag Zo in 13.61 we take the zigzag in [Z] 
with the highest index. Note that while performing repairing move 3 
we have first deleted from the matrices Iz, Pz and Qz all rows which 
corresponded with intersection points on one of the zigzags in the sequence 
(Zi, . . . , Zis)- Subsequently we have inserted rows for intersection points 
of a zigzag in the new sequence (Z(, . . . , Z' 2a ) with a zigzag which also 
intersects Z' = Zq. The previously applied repairing moves 2 had already 
removed all intersection points of Zo with the zigzags in [Z] which were 
not in (Zi, . . . , Zgs). So after applying repairing moves 3 two zigzags in 
[Z] U (-[Z]) do not intersect. 

4.11. After the above merging and repairing moves we have produced 
a very good pattern of zigzags and now return to 14.61 for the next 
iteration. 

5 From Bz, Iz, Pz, Qz to Kg and back 

5.1. The conversion works for every good pattern Z of, say p, zigzags. 
Condition 12.121 9 is not needed here. The rows of Pz and Qz must be 
taken modulo the row space of Bz- This is achieved by multiplying Pz 
and Qz from the right by a p x (p — 2)-matrix Az with entries in Z, such 
that rank(Az) =p-2 and Bz A z = 0. 

The rows of Pz Az represent points in V/l^Bz- We denote the set of 
these points by 03, because these are in fact the black nodes in the dimer 
model. In the zigzag pattern these are the +-cells. Similarly, we denote 
the set of rows of Qz Az by 2U. These are the white nodes in the dimer 
model and the — cells in the zigzag pattern. 
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5.2. Definition. (cf. 5 Definition 8.2, 6 Definition 1) The generalized 
Kasteleyn matrix Kz(z, u) of a good pattern of zigzags Z — (Zi, . . . , Z p ) 
is defined as follows. The rows of Kz(z,u) correspond 1 : 1 with the 
elements of 25 and the columns correspond 1 : 1 with the elements of 23. 
The entries of Kz(z, u) are polynomials in two sets of variables z and u. 
The variables in u = {iii, . . . , u p } correspond 1 : 1 with the zigzags in Z, 
and, hence, with the columns of Bz, Iz, Pz and Qz- The variables in 
z = {z e } correspond 1 : 1 with the intersection points e of Z and, hence, 
with the rows of Iz, Pz and Qz- 

For an intersection point e we denote by b(e) the element of 25 which 
"is" the e-th row of PzAz- Similarly, w(e) £ 2H "is" the e-th row of 
Qz Az- Finally, i(e) and j(e) are such that Iz(e,i(e)) — Iz(e,j(e)) = 1. 

Finally, for b 6 23 and w £ 2D we define: 



the (b,-w)-entry of Kz(z,u) 



e: b(e) — b, w(e) — w 



(11) 



5.3. Example. For the pattern of zigzags in Figure [2] the generalized 
Kasteleyn matrix is: 



K z (z,u) = 



21U1M5 Z2U1U2 + Z3U4U5 Z4U2U4 

25M1U3 Z6^lU6 + Z7U3U4 Z8M4U6 
Z9U3U5 Z10M5U6 + Z11U2W3 2l2«2«6 



5.4. Definition. 

Kasteleyn matrix 



(cf. [B] Definition 3) The complementary generalized 
S-z{ z i u ) °f a good pattern Z of p zigzags is 



vj(z,u) 



ui 



5.5. Remark. The information in Kz(z, u) is in fact equivalent with 
that in Bz, Iz, Pz- By Theorem 9.3 and §3.7 in [5] the columns of 
Bz are the primitive vectors along the sides of the Newton polygon of 
det Kz(z, u) w.r.t. ui, . . . , u p . One recovers Pz and Iz as follows. Let V 
and Q, respectively, be the sets of exponent vectors of the monomials in 
ui , • ■ • , Up which appear in the first columns of the matrices 

K z (z,u) (K^(z,u)K 2; (z,u))' 1 resp. (K^(z,u)K z (z,u))" +1 

for n £ Z>o. Translations on IF by vectors in T?Bz preserve V and Q. 
Now take a 'fundamental domain' V* C V for the Z 2 _B^-action on V ■ For 
each vector a in V* let Q a be the set of vectors j3 in Q such that a — (3 
has precisely two non-zero entries and these are both 1. The rows of Pz 
resp. Iz are a resp. a — f3 with a £ V* and f3 £ Q a . 



6 The principal ^-determinant 

6.1. The principal ./La-determinant for a good pattern of zigzags. 

Let Z — (Zi, . . . , Z p ) be a good pattern of zigzags and let Az denote the 
set of rows of the matrix Az fsee !5.1[) . Let t be the homomorphism 

1 : Z[z e I e intersection point in Z] — > Z , t(z e ) = |^( e ) A ^j( e ) | • 
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The principal Az -determinant is defined in [2] and Theorem 3 in [6] states 
that it is equal to 

t(detK^(z,u)) . (12) 

The relation 

Iz = Pz - Q z 
is precisely the one required in Condition 2 of |B]. 

6.2. Concluding remark about the principal ^-determinant. 

For the principal ,4-determinant of the set A = {ai, . . . ,ajv} in the In- 
troduction we may have to make some slight adaptations to formula (112[) . 
which reverse in a sense the transformation from Ba to Bz in l4.ll 

In (|12p the variables in u = (ui,...,u p ) correspond 1 : 1 with the 
columns of Bz- Take a new set of variables v = (v\, . . . ,vm) which 
correspond 1 : 1 with the columns of Ba- Recall that Bz is obtained 
from Ba by permuting and splitting up columns as in Q. Then, to 
reverse the transformation from Ba to Bz one must set m = dk Vk if the 
i-th column of Bz comes from the fe-th column of Ba', here dk is the g.c.d. 
of the two entries in the fe-th column of Ba- 
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